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ABSTRACT 

We present the broad-band Spectral Energy Distributions (SEDs) of the largest 
available highly (72%) complete spectroscopic sample of mid-infrared (MIR) se- 
lected galaxies and AGN at intermediate redshift. The sample contains 203 
extragalactic sources from the 15-/im survey in the ELAIS-SWIRE field SI, all 
with measured spectroscopic redshift. Most of these sources have full multi- 
wavelength coverage from the far-UV (GALEX) to the far- infrared (Spitzer) and 
lie in the redshift range 0.1 < z < 1.3. Due to its size, this sample allows us for 
the first time to characterise the spectral properties of the sources responsible 
for the strong evolution observed in the MIR. Based on SED-fitting technique 
we have classified the MIR sources, identifying AGN signatures in about 50% of 
them. This fraction is significantly higher than that derived from optical spec- 
troscopy (~29%) and is due in particular to the identification of AGN activity 
in objects spectroscopically classified as galaxies. This might be partially due to 
the fact that the spectroscopic classification can be somewhat unreliable because 
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of host galaxy dilution in the optical. It is hkely that in most of our objects, 
the AGN is either obscured or of low-luminosity, and thus it does not dominate 
the energetic output at any wavelength, except in the MIR, showing up just in 
the range where the host galaxy SED has a minimum. The fraction of AGN 
strongly depends on the flux density, with that derived through the SED-fltting 
being about 20% at 5*15^^ ~ 0.5 — 1 mJy and gradually increasing up to 100% 
at Si^^rn > 10 mJy, while that obtained from optical spectroscopy never being 
>30%, even at the higher flux densities. The results of this work will be very 
useful for updating all the models aimed at interpreting the deep infrared survey 
data and, in particular, for constraining the nature and the role of dust-obscured 
systems in the intermediate/high- redshift Universe. 

Subject headings: infrared: galaxies - galaxies: evolution - galaxies: Seyfert - 
galaxies: starburst - cosmology: observations. 

1. Introduction 

Understanding the overall Spectral Energy Distribution (SED) of sources responsible for the 
observed Cosmic Infrared Background (CIRB) is a crucial tool for unveiling the nature and 
evolution of galaxies and AGN in the infrared (IR) and for obtaining a complete picture 
of the history of star-formation and obscured AGN activity in the Universe. The mid- and 
far-infrared (MIR and FIR) regions of the electromagnetic spectrum efficiently probe the 
population of actively star-forming galaxies and dust-obscured AGN. 

The InfraRed Astronomical Satellite IRAS has sampled the local Universe {z < 0.2) in 
the MIR/FIR band, discovering Ultra Luminous Infrared Galaxies (ULIGs: Ljji > lO^^L©; 
Sanders & Mirabel 1996, Lonsdale, Farrah & Smith 2006) and flrst showing some evidence 
for their strong evolution (i.e. Hacking, Houck and Condon 1987, Lonsdale & Hacking 1989; 
Saunders et al. 1990). 

Extragalactic source counts from different surveys over a wide flux range obtained with the 
/.S'O CAM instrument (Cesarsky et al. 1996) on board of the Infrared Space Observatory {ISO; 
Kessler et al. 1996) indicate that also IR sources less extreme than ULIGs (i.e. Luminous 
Infrared Galaxies; LIGs: Ljr > IO^^Lq) have evolved rapidly from z ~ up to z ~ 1.2, 
significantly faster than deduced from optical surveys (i.e. Elbaz et al. 1999; Gruppioni et 
al. 2002; Pozzi et al. 2004). 

The Spitzer Space Telescope (Werner et al. 2004) is now providing new insights into the 
IR population contributing to the CIRB, in particular in the MIPS 24-/xm band, where the 
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extragalactic source counts (Marleau et al. 2004; Papovich et al. 2004) have confirmed the 
existence of the rapidly evolving dust-obscured population discovered by ISOCAM, detecting 
the higher-^ analogues (1 < 2; < 3) of the /S'OC^M galaxies. In particular, the deeper MIPS 
surveys are now able to resolve about 70% of the CIRB at 24 fira, with the bulk of this 
background originating in discrete sources ai z^2 (e.g. Papovich et al. 2004). 

These results have stimulated the development of several evolutionary models for IR 
galaxies (i.e. Rowan-Robinson 2001; Franceschini et al. 2001; Chary & Elbaz 2001; Xu et al. 
2003; Lagache, Dole & Puget 2003; Lagache et al. 2004; Pozzi et al., 2004; Gruppioni et al. 
2005; Pearson 2005), aimed at reproducing the MIR source counts and the observed rcdshift 
distributions. All these models are more or less able to reproduce the observed ISOCAM 
15 iim source counts, but none of them can provide an acceptable fit to the MIPS 24-/xm 
counts without any "ad-hoc" changes. In particular, they tend to predict the characteristic 
turn-over of the differential 24- //m number counts at fiuxes higher than observed. Even worse 
is the situation regarding the observed redshift distributions, both for ISOCAM and MIPS 
sources. Le Floc'h et al. (2005) and Perez-Gonzalez et al. (2005), by studying the luminosity 
evolution of a sample of MIPS 24-/im sources in the Ghandra Deep Field South (GDFS) up 
to z = 1 and in the range 1 < z < 3 respectively, find that all the models considered in 
these works (Lagache et al. 2004, Ghary et al. 2004, Gruppioni et al. 2005 and Pearson et 
al. 2005) fail in reproducing the differential counts at 24 /xm in different redshift shces and 
the total ^-distribution. We must note that most of the redshifts considered in these works 
are photometric and the analysed fields are relatively small, and can therefore be affected 
by cosmic variance. However, there seems to be a significant disagreement between all the 
evolutionary models existing in literature and the recent Spitzer data. Two main weaknesses 
are common to all these models: 

1) They are based on very uncertain assumptions about the shape of the SED of MIR 
sources. In particular, they generally assume an invariant PAH emission and extrapo- 
late few local SEDs to high redshifts and luminosities (L), either by keeping them fixed 
or artificially increasing their activity with luminosity (e.g. Ghary & Elbaz 2001). 

2) They severely underestimate, or even completely neglect (e.g. Lagache et al. 2003; 
Ghary & Elbaz 2001), the AGN contribution in MIR sources (see Brand et al. 2006; 
Gruppioni et al. 2006). In fact, without MIR spectroscopic data or, if not available, 
without a complete SED characterisation, it is very difficult to disentangle star-forming 
galaxies from AGN (Genzel et al. 1998, Sajina et al. 2007). This is particularly true in 
cases where spectroscopic classification based on optical line diagnostics is impossible 
due to the absence of some of the required emission lines in the observed spectral 
range, or to the significant differential line extinction produced by dust, or to the lack 
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of spectroscopic data and deep X-ray coverage. Since different evolutionary models are 
applied to the two populations, the misidentification of AGN and star- forming galaxies 
introduces large uncertainties in the predicted source counts and backgrounds.. 

The crucial wavelength range (3-200 /im) covered by Spitzer with unprecedented sensi- 
tivities gives the opportunity to characterise for the first time the MlR/FlR spectral prop- 
erties of large numbers of sources over a significant extent in redshift and to study their 
evolution with z and/or L, testing the model assumptions. In particular, thanks also to the 
extensive multi-wavelength coverage available in several areas of the sky, it is now possible to 
study the broad-band (from UV to FIR) spectral properties of the same sources responsible 
for the observed evolution in the MIR, thus constructing observational templates libraries 
over a large range of wavelengths. 

In the Southern field of the European Large Area ISO Survey (ELAIS; Oliver et al. 
2000; Rowan-Robinson et al. 2004), SI, an extensive multi-wavelength follow-up campaign 
has been performed during the past years. This field, initially selected for ISO observations 
at 7, 15 and 90 /xm, is now one of the six areas covered by the Spitzer Wide- Area InfraRed 
Extragalactic Survey (SWIRE, Lonsdale et al. 2003) in all the IRAC and MIPS bands. 
In addition to ISO and Spitzer observations, the whole SI field (~4 deg^) has been deeply 
surveyed in the radio {ATCA; Gruppioni et al. 1999; Ciliegi et al., in preparation), in the 
optical {R band, ESO/3.6-m; La Franca et al. 2004) and in the near- and far-UV [Galaxy 
Evolution Explorer (GALEX); Burgarella et al. 2005]. The central square degree of SI is 
covered by NIR (ESO/SOFI; Matute et al. in preparation) and optical {ESO/WFI; Berta et 
al. 2006) data, while ~0.6 deg^ have been observed in the X-rays {Beppo-SAX; Alexander 
et al. 2001 and XMM-Newton; Puccetti et al. 2006). Optical spectroscopic informations 
are available for ~200 15-/im ISOCAM sources {ESO/3.6-m; La Franca et al. 2004) and for 
most of the XMM and some of the SWlRE-24/im sources having optical counterpart brighter 
than R ^ 24 (with ESO/VIMOS; Feruglio et al., 2008; La Franca et al., in preparation). 

Due to its extensive multi-wavelength coverage and its large area extent, SI is well 
suited to study in detail the broad-band SEDs of a statistically significant sample of infrared 
galaxies and AGN at intermediate redshifts (z < 1.5). To this purpose, we have selected the 
203 extragalactic sources detected by ISOCAM at 15 /im {S > 0.5 mJy) in the SI field (Lari 
et al. 2001) with R < 23, available spectroscopic redshift and classification (La Franca et al. 
2004; La Franca et al. 2007). The same sample has previously been used to derive the first 15- 
//m luminosity function of galaxies (Pozzi et al. 2004) and AGN (Matute et al. 2002, 2006). 
The treatment of infrared/optical incompleteness of this sample for statistical purposes is 
described in detail by Gruppioni et al. (2002) and La Franca et al. (2004). The sample 
discussed here can be considered as the largest available spectroscopic sample of IR galaxies 
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and AGN at intermediate redshift, with such a high level of completeness (72%), allowing 
us for the first time to characterise the spectral properties of the sources responsible for the 
strong evolution observed in the MIR. Although large spectroscopic samples are available 
in fields observed by Spitzer, like the GOODS-CDFS (see Vanzella et al. 2005; 2006), the 
targets are not selected to be MIR sources, but optical/NIR objects responding to given 
colour criteria and likely to be at high redshifts. In the very near future large spectroscopic 
samples of MIR (24-/im) selected sources, either locally (i.e. in the FLS by Marleau et al. 
2007) or at high redshift (i.e. in the COSMOS area) will be made available to the community, 
allowing studies similar to the one presented here, but at present our sample is unique to this 
purpose. Indeed most of our objects are LIGs, with MIR fluxes in the crucial range between 
the IRAS surveys (200 mJy) and the Deep /50C^M Surveys (0.1 mJy; Elbaz et al. 1999), 
where MIR source counts start diverging from no-evolutionary expectations. Most of these 
objects (except for type 1 AGN) are in the redshift range 0.1 < z < 1.0. We make use of 
all the available data, from FUV to FIR, to derive the SEDs of these sources and construct 
the first observational library of templates for MIR galaxies and AGN at intermediate z. To 
interpret the observed SEDs we perform a fit with several local template SEDs, representative 
of different classes of IR galaxies and AGN (PoUetta et al. 2007), comparing the resulting 
SED classification with the spectroscopic one. 

The present paper is structured as follows. In Section [2] we describe the reference sample. 
In Section [3] we present the multi-wavelength data-set and the observed SEDs, discussing 
the SED-fitting procedure and results. In Section H] we compare the original spectroscopic 
classification with the SED classification. In Section [S] we present the revised 15-/im counts 
for AGN and galaxies, we discuss the results and present our conclusions. 

Throughout the paper we adopt a cosmological model with Qrn = 0.3, = 0.7 and 
Hq = 70 km s"^ Mpc^^. The magnitudes discussed here are in the Vega system. 

2. The Sample 

Our reference sample is the Lari et al. (2001) catalogue of 462 sources selected at 15 fim 
(LW3 band) in the ELAIS Southern field SI (centred at a(2000) = OO'^ 34™ 44.4^ 5(2000) = 
-43° 28' 12", and covering about 2° x 2°). This sample, complete at the 5a level, is the only 
ISO CAM sample covering the whole flux density range 0.5 - 150 mJy, thus linking IRAS to 
the Deep ISOCAM Surveys (Elbaz et al. 1999). The source counts at 15 fim obtained from 
that catalogue sample the flux density region where observed counts start diverging from no 
evolution models, as discussed by Gruppioni et al. (2002). La Franca et al. (2004) presented 
R-band data for a highly reliable sub-sample of 406 out of the total 462 15-/im sources of the 
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Lari et al. (2001) catalogue. The R band data were obtained by an ESO imaging campaign 
with the DFOS'C instrument mounted on the 1.5-m ESO/Danish telescope at La Silla (Chile), 
providing a reliable optical counterpart down to i? ~ 23 for 317 of them, thus reaching 95% 
completeness level. Spectroscopic observations of the optical counterparts of the ISO CAM 
SI sources were carried out at the AAT/2dF, ESO/Danish 1.5-m, 3.6-m and NTT telescopes 
(La Franca et al. 2004), providing a secure spectroscopic identification for 290 ISOCAM 
sources (72% of the whole highly rehable ELAIS-Sl sample). La Franca et al. (2004) have 
classified 199 of the spectroscopically identified sources as extragalactic objects (25 type 1 
AGN, 23 type 2 AGN, 9 hners, 32 starburst galaxies, 100 Ha emitter galaxies, 3 early type 
galaxies, 7 unclassified, but with measured redshift), the remainder 91 as stars. Details about 
the identification completeness of the extragalactic sample for statistical uses are given in La 
Franca et al. (2004). These 199 extragalactic objects with measured redshift are the same 
objects previously used for the statistical analysis of the evolution of galaxies and AGN 
by Pozzi ct al. (2004) and Matute et al. (2002, 2006) respectively. With respect to the 
La Franca et al. (2004) catalogue, two more 15 /xm sources were spectroscopically identified 
through VIMOS-VLT observations (La Franca et al. 2007): ELAISG15_J003317-431706 is a 
R = 24.3 galaxy showing [Oil] emission at redshift 0.689, while ELAISG15_J003447-432447 
is an AGN2 at redshift 1.076. Moreover, we derived redshifts for two additional 15 /xm 
sources: for ELAISC15_J003915-430426 we found a z value of 0.013 in the NED database, 
while for ELAISC15_J003545— 431833 we were able to measure z through a more accurate 
reduction of the spectrum. Three sources with previously poor quality spectra changed 
their spectroscopic classification after having been re-observed with VIMOS at ESO-VLT 
(La Franca et al. 2007): ELAISC15_J003330— 431553, which was wrongly classified as a 
starburst galaxy z — 0.473, showed broad GUI and Mgll emission dX z — 2.170 typical 
of AGNl activity; ELAISG15_J003603-433155, which was classified as AGN2, showed a 
broad Mgll emission typical of AGNl activity; and ELAISG15_J003622-432826, which was 
classified as a starburst galaxy, showed a clear [0III]/H/3 ratio typical of AGN2 activity. 

In summary, we have considered 203 extragalactic sources selected at 15 //m with R < 
23, spectroscopic redshift and multi-band photometry. 

3. The SEDs 

3.1. The Multi- Wavelength Data-Set 

The objects we focus on cover the redshift range 0.01 - 3.09 (0.01 - 1.30 if excluding type 
1 AGN) and more than two orders of magnitude in 15-/im flux (0.5 - 60 mJy). The multi- 
band SEDs for these sources have been constructed by looking for counterparts in all the 
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multi- wavelength catalogues from far-UV to FIR available in the ELAIS-Sl field. Given for 
granted the optical R-band/spectroscopic association performed by La Franca et al. (2004), 
we have searched for associations in the GALEX Deep Imaging Survey (DIS), ESIS, J,Ks 
and SWIRE catalogues. The cross-correlation radii between the 15-/im or optical positions 
and the different matching catalogues positions have been chosen equal to the value above 
which spurious associations start dominating (i.e. Hook et al. 1998). 

The ultraviolet Galaxy Evolution Explorer {GALEX, Martin et al. 2005) Deep Imaging 
Survey (DIS) is observing ~100 deg^ in 12 different areas of the sky, including ELAIS-Sl 
(Burgarella et al. 2005). Twelve ELAIS-Sl G^dLi^X tiles are already available to the public as 
part of the second data release (GR2). The exposure times for these tiles vary from ~3000 
to ~50000 sec. The photometric catalogues in the far-UV (FUV; 1530 A) and near-UV 
(NUV; 2310 A) have been cross-correlated with the optical positions of the ISOGAM source 
counterparts, with a matching radius of 3 arcsec. Then all the possible UV counterparts 
have been visually inspected on the GALEX images, in order to reject spurious associations. 
We found 176 likely associations with the NUV catalogue, 169 of which are also in the FUV 
one. 

The ESO-Spitzer Imaging Extragalactic Survey (ESIS; Berta et al. 2006) has covered the 
central 1.5 deg^ of ELAIS-Sl with WFI/2.2-m deep observations in the B, V and R bands, 
reaching 95% completeness a.t B,V ^ 25 and R ~24.5. The cross-correlation between the 
original optical positions and the ESIS positions (within 2 arcsec) have produced a good 
match for 90 sources (those falling in the area covered by ESIS), thus providing additional 
optical {B, V) magnitudes for about 45% of the 203 /S'OC^Msources. In the area not covered 
by ESIS, R magnitudes from the La Franca et al. (2004) catalogue and bj magnitudes from 
APM survey (~90% complete to bj = 20.5; Maddox et al. 1990) have been used. 

The central square degree of the ELAIS-Sl area (Sl-5) was the target of NIR J and Ks 
imaging with NTT/SOFI, reaching J ~21 and Ks ~20 (Vega)0 (Matute et al., in prepa- 
ration). The optical positions of the 15-/im source counterparts have been cross-correlated 
with the NIR catalogue ones, using a matching radius of 2 arcsec. Thirty unique NIR coun- 
terparts have been found in the central square degree of ELAIS-Sl, while additional NIR 
data from the Two Micron All Sky Survey (2MASS; Jarrett et al. 2000) are available over 
the whole ISO area, providing 110 matches (10 of which are in common with the Sl-i^s 
ones) . In total we found s-band magnitudes for 130 and J-band magnitudes for 30 out of 
the 203 15-/im sources. 

The Spitzer/SWIKE observations in ELAIS-Sl cover the whole ISO region and a total 
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sky area of ~7 deg^ in all the IRAC and M/PS* bands, reaching 5a sensitivities of 4.1, 8.5, 48.2, 
53.0, 256.0, 26x10^ and 166x10^ /xJy in the 3.6, 4.5, 5.8, 8.0, 24, 70 and 160 jjm channels 
respectively (Lonsdale et al. 2004). The SWIRE data in ELAIS-Sl have been released to the 
community in autumn 2005 (Data Release 3, DR3) through the query building GATOR at the 
NASA/Infrared Science Archive (http : // irsa. ipac . caltech. edu/ applications/Gat or/ ). 
However, to the purposes of this work, we had access to the SWIRE working catalogues, 
which contain all of the sources in the public catalogue, but reach deeper flux densities. De- 
tails about the SWIRE data reduction, that was carried out by the Spitzer Science Center 
and SWIRE team, can be found in the Data Release paper (Surace et al. 20of|). For the 
IRAC fluxes and positions we used the band-merged working catalogue, consisting of IRAC 
and MIPS-24 fim fluxes associated with each other. While the public catalogue is a 3.6-/im 
based one, requiring detections at 3.6 and 4.5 fim above specific SNR thresholds (10 at 3.6 
fim and 5 at 4.5 /im), the working catalogue contains all detections, at any signal-to-noise 
ratios and even if a source is detected only in one IRAC hand. For the 24-, 70- and 160-/im 
fluxes and positions, we used the single-band only catalogues, which contain all the sources 
detected only in the MIPS waveband in question. The IRAC positional accuracy is about 
0.2"-0.4" in all bands, while the MIPS ones are ~1", 4" and 8" at 24, 70 and 160 /im re- 
spectively. Fixed aperture photometry is used for point-like objects, corrected for aperture 
losses (we have considered the SWIRE "aperture2" (1.9 arcsec in radius) for both IRAC 
and MIPS), while "ad-hoc" aperture photometry (i.e. as large as needed to recover the 
whole flux of very extended sources) has been performed by us directly on images for the 
extended sources, for which even the SExtractor's Kron fluxes (those suggested for extended 
sources) are found to underestimate the total flux densities (up to a factor of ~2 at 24 
fxm). For the match with the /i?j4 C band-merged catalogue we have considered the optical 
positions, searching for counterparts within a radius of 3 arcsec, while for match with the 
MIPS catalogues we have considered the ISOCAM positions with searching radii of 5, 10 
and 20 arcsec at 24, 70 and 160 fim respectively. All but one of the 203 15-/im sources are 
detected in at least two /i?74C bands (most in all 4 bands), while 200 are detected at 24 
/im, 110 at 70 fim and 59 at 160 fim. The matches with 24-, 70- and 160-/im sources have 
been visually inspected to be sure they were not mismatches. From the visual inspection of 
the MIPS images we have found clear detections for 43 and 45 additional sources at 70- and 
160-/xm respectively. However, we have chosen not to derive flux densities for these sources, 
considering in this work only those included in the SWIRE catalogues. 

In Table 1 we report the multi-wavelength informations available for all the 203 15-/xm 
sources. The source name from the Lari et al. (2001) catalogue is reported in the first 



also http : //data . spitzer . caltech . edu/popular/swire/20050603_enhcinced_vl/Docuinents/SWIRE2 jloc_083105 . pd: 



- 9 - 



column, while the flux density (in fiJj) in all bands, from FUV to FIR, the redshift and 
the total infrared luminosity (obtained by integrating the best-fltting SED in the range 8 - 
1000 fim) are reported in the following columns. In the last two columns the spectroscopic 
and SED- fitting classifications (see Section H]) are presented. In case of no detection we have 
quoted a 3a upper limit, while — means "no data coverage". 

3.2. The Template Library 

We have made use of all the available data, from FUV to FIR, to derive the SEDs of 
these 203 sources and construct the first highly complete observational sample of SEDs for 
MIR galaxies and AGN at intermediate z to be used for statistical purposes (i.e. luminosity 
function). To interpret and classify the observed SEDs, we have performed a fit with several 
local template SEDs, representative of different classes of IR galaxies and AGN (including 3 
Ellipticals of different ages, 1 lenticular, 7 Spirals, 3 Starbursts, 3 QSOs, 1 reddened QSO, 
Seyfert 1, 1.8 and 2 and 2 ULIRGs, containing both a starburst and an AGN component, 
in the wavelength range between 0.1 and 1000 fim (Polletta et al. 2007), appropriately 
redshifted to the rest frame. The full library of galaxy and AGN templates (excluding the 
three ellipticals since none of our sources could be fitted by any of them) is shown in Fig. [TJ 
The elliptical, spiral and starburst templates were generated with the GRASIL code (Silva 
et al. 1998). The spiral templates range from early to late types (SO-Sdm), the starburst 
templates correspond to the SEDs of NGC 6090, M 82 and Arp 220. Templates of moderately 
luminous AGN, representing Seyfert 1.8 and Seyfert 2 galaxies, were obtained by combining 
models and data from a large sample of Seyfert galaxies. The AGN templates include 
three templates representing optically-selected QSOs with different values of IR/optical flux 
ratios (QSO, QSO_high, and QSOJow) and one reddened QSO (red_QSO). The composite 
(AGN-I-SB) templates are empirical templates created to flt the SEDs of the heavily obscured 
BAL QSO Mrk 231 (Berta 2005) and the Seyfert 2 galaxy IRAS 19254-7245 South (Berta 
et al. 2003). These objects contain a powerful starburst component, mainly responsible 
for their FIR emission, and an AGN component that contributes to the MIR (Farrah et al. 
2003). With respect to existing template libraries derived from empirical SEDs (Coleman et 
al. 1980) or from models (Bruzual 2003; Fioc et al. 1997; Silva et al. 1998; Devriendt 1999), 
this library has a broader wavelength coverage and variety of spectral types. Examples of 
application of this library to various types of SWIRE sources can be found in Lonsdale et 
al. (2004), Franceschini et al. (2005), Hatziminaoglou et al. (2005), Jarrett et al. (2006), 
Polletta et al. (2006), Weedman et al. (2006), Tajer et al. (2007), Polletta et al. (2007). 

To avoid biasing the results in favour of a small subset of templates, we adopt the fuU 
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library of Polletta et al. (2007) to fit the SEDs. Although some templates are obtained 
from modelling the SEDs of Ultra-Luminous objects (ULIGs, Lui > IO^^Lq: e.g. Mrk 231, 
119254 and Arp 220), they are kept in the library since they well represent a large fraction 
of the Spitzer population at several flux and luminosity levels (see e.g. Franceschini et al. 
2005, Polletta et al. 2006, 2007, Alonso-Herrero et al. 2006). 

3.3. The Fitting Method 

The SEDs of our sources are fitted using an automated fitting routine contained in the 
Le PHARE software (by Arnouts & Ilbert, available at |http : //www . lam . oamp . f r/people/| 
arnouts/LE_PHARE.html), offering the possibility of using any template library and filters 
set. Le PHARE (PHotometric Analysis for Redshift Estimations; Ilbert et al. 2006) is a 
publicly available set of fortran programs aiming at computing photometric redshifts through 
the best-fitting SED analysis. The program is based on a simple fitting method between 
the expected magnitudes and an observed photometric catalogue. 

Since we already know the spectroscopic redshift of our galaxies, we have fixed z and have 
used Le PHARE only to obtain the best-fitting SED through a comparison with the template 
SEDs. No extinction was added to the templates, since some of them are already intrinsically 
extinct. Note that by using a spectroscopically complete sample we are able to avoid all 
the uncertainties due to the photometric redshift measurements, which greatly increase the 
degeneracy in the SED-fitting procedure. Moreover, for the first time we can compare the 
results from optical spectroscopy with the broad-band SED shape of a statistically significant 
sample of MIR galaxies and AGN. 

We make the basic assumption that the SED shapes seen at low redshifts can also well 
represent the higher redshift objects. At the end of our procedure we would be able to 
verify this assumption. The template library used to fit our data contains a finite number 
of SEDs (21), representative of given classes of local infrared objects, which do not vary 
with continuity from one class to another (there are large gaps in the parameter space). 
Therefore, the quality of the fit depends not only on the photometric errors, but also on 
the template SED uncertainties. For this reason, in our fitting procedure, in addition to 
the photometric errors on data, we need to take into account also the uncertainties due to 
the template SEDs "discretization" . To estimate and consider the total uncertainty on both 
the photometric data and the template SEDs, we have proceeded as follows. First, we have 
increased all the formal catalogue errors (probably too small) by a reasonable given amount, 
going from 5% in the UV, and optical bands up to ~15% in the MIPS bands. We have run 
Le Phare on all the sources, obtaining a first "estimate" of the best-fitting SED classification. 
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We have then computed the distributions of the {Sobject - S template) band / {cr object) band values in 
each of the considered photometric band (where Sobject and aobject are the flux density and 
the relative error of the source, and Stempiate the flux of the template in the considered band), 
iteratively increasing the photometric errors until we have obtained a Gaussian distribution 
with cr ~ 1. This corresponds to reduced distributions peaked around 1 (as should be 
in case of good fit). The enlarged relative errors {cr object/ Sobject) obtained in this way go, on 
average, from ~12-15% in the optical/NIR//i?.4C bands up to ~20% in the MIPS hands. 
With the new (significantly larger) photometric uncertainties, which now take into account 
also the additional uncertainties induced by the template SED discretization, we have run 
Le PHARE on our sources for the second time, obtaining what we have taken as the final 
SED-fitting results. 



3.4. The SED-Fitting Results 

In figure [2] examples of data overlying the best-fitting template SEDs are shown for 
most of the templates reproducing our data (the original library contains 21 SEDs, but none 
of our sources could be fitted by any of the 3 elliptical galaxy templates). We find that 83 
sources (41%) are reproduced by a galaxy template SED (4 SO, 3 Sa, 16 Sb, 27 Sc, 18 Sd, 15 
Sdm), 13 (6.5%) with a starburst template SED (1 Arp 220, 5 M 82, 7 NGC 6090), 72 (35%) 
with a type 2 AGN template SED (52 Seyfert 2, 15 Seyfert 1.8, 5 red QSO), 11 (5.5%) with 
a ULIG composite SED (9 Markarian 231, 2 119254) and 24 (12%) with a type 1 AGN SED 
(3 Seyfert 1, 6 QSO, 13 high IR/optical flux ratio QSO, 2 low IR/optical flux ratio QSO). 
Therefore, with our broad-band SED-fitting method we find AGN activity in 52.5% of our 
MIR selected sample. 

We note that the distributions of the final logm^S object/ Stempiate) values in some bands 
and for some SED classes are centred on values slightly different from 0. These offsets, 
which are different in the different bands and for the different best-fitting SED classes (i.e. 
for the QSO class the larger offsets are observed at 70 and 160 /xm, while for the galaxy 
class the larger offsets are observed in the B, J, K, 15-/im and 24-/im bands; see Fig. [3]), are 
interpreted as the average values we should add to the template SEDs to better reproduce 
the data ensemble of MIR sources. In Fig. [3] we show the rest-frame SEDs (dots) of 15-/im 
sources belonging to the different "broad" SED classes (AGNl, ULIG, AGN2, Starburst and 
Spiral galaxy), compared to the most frequent best-fit template SED of that class (i.e. "QSO" 
for AGNl, "Seyfert 2" for AGN2, "Mrk 231" for ULIG, "NGC 6090" for Starbust, "Sc" for 
Spiral) , all normalised to the s-band flux density. The white triangles represent the average 
values that the templates should have in the different rest-frame bands to better reproduce 
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the bulk of our sources, while the grey-shaded areas represent their relative uncertainty 
region corresponding to ±la of the expected values. From this check we can conclude 
that, on average, the template SEDs of the PoUetta et al. (2007) library modelled on local 
galaxies and AGN are able to reproduce also most of the MIR selected sources at intermediate 
redshifts, though some of the local templates should be slightly modified in some wave-bands 
to get a better agreement with the bulk of the observed data. 

In order to estimate the reliability of the SED classification, we compared the classi- 
fication and obtained by the two best solutions of each fit: the primary (corresponding 
to the minimum value, xlest) secondary solutions (xLd)- ^^^^ P"" 

mary solutions of our photometric analysis are quite stable within uncertainties. In fact, for 
the majority of our sources (92%) the xlest values are within the 90% probability threshold 
(^90% ~ 21.07 for the case of 14 degrees of freedom). In most cases the secondary solution 
belongs to the same "broad" SED class of the primary (i.e. AGN or galaxy, containing all the 
type I's, type 2's and ULIG templates and all the galaxy and starburst templates, respec- 
tively) and only for 23 objects (11%) it belongs to a different class, with a xlnd lower than 
the xlo% threshold. In all the other cases, if the primary and secondary solutions correspond 
to different SED classes, the secondary solutions have a low probability of being acceptable 
ixlnd > X9o%)- The 16 sources (8%) with a statistically bad primary solution {xlest > xloyj 
have, in a few cases, some evident photometric errors in one or more bands (mainly in the 
J, K' ones) . In all the remaining cases, the bad fit is due to an observed SED flatter than the 
best-fitting template in the MIR domain, thus implying the need of a template SED (missing 
in our library) intermediate between the power-law and the Seyfert 2/1.8 ones, where the 
AGN is not dominant as in the type 1 template case (power-law), but its contribution is 
higher than in the Seyfert 2 template case. 



4. Spectroscopic versus SED Classification 

We have compared the results of our SED-fitting procedure with the results obtained 
from the optical spectroscopy. In Table [2] we show for each spectroscopic class the number 
of sources which are best fitted by the various SED templates and vice-versa. The reported 
percentages are the fractions of sources of a given spectroscopic classification fitted by the 
various template SEDs. The spectroscopic type 1 (broad hne) AGN are all fitted by AGN 
template SEDs, and therefore the agreement between the two classifications for type 1 AGN 
is very good (the only three objects fitted by a type 2 AGN template are indeed fitted by 
a Seyfert 1.8 template). However, the agreement is not equally good for the spectroscopic 
type 2 AGN (+ LINERs) and galaxies: 56% of the AGN2+LINER sample is fitted by a type 
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2 AGN template and another 19% is fitted by type 1 AGN or composite (AGN+starburst) 
templates, but ~25% do show a galaxy SED. Also for the majority (63%) of the spectroscopic 
normal galaxies the two classifications agree (another 3% is fitted by a starburst template), 
but there is a significant fraction (34%) fitted by AGN templates (mostly type 2's). For the 
spectroscopic starbursts the fraction of sources classified as AGN by the SED fitting is even 
higher, 53%. The small number of spectroscopically unclassified sources is fitted for the 40% 
by AGN and for the 60% by galaxy templates. Therefore, we can conclude that the SED- 
fitting technique apphed to our MIR selected sample is able to identify AGN activity in ~40% 
of sources spectroscopically unclassified, in ~34% of sources spectroscopically classified as 
galaxies, in ~53% of spectroscopic starbursts, in 75% of type 2 AGN or LINERs, and in 
100% of type 1 AGN. The main result of this comparison is that, although for many sources 
the spectroscopic classification is in agreement with the classification resulting from the 
SED-fitting method, the broad-band SED analysis is able to find AGN activity in a higher 
fraction of MIR sources (~53%) than the optical line diagnostic techniques (~29%). On the 
other hand, a small number (14%) of the spectroscopically classified AGN are classified as 
galaxies by the SED fitting. 

For comparison, we cite the recent results of Polletta et al. (2007), who, using a sample 
of X-ray selected AGN with available spectroscopic classification and a similar fitting method 
as applied to our sample, find that 82% of the X-ray selected AGNs are classified as AGN by 
the SED-based classification method, compared to 78% by the spectroscopic classification. 
However, the SED and the spectroscopic classification agree only in 53% of the cases: the 
SED method and the spectroscopic classification for X-ray sources are consistent for all type 
1 AGN, but for only 29% of type 2 AGN and 33% of AGN with galaxy spectra. 

The more evident change between the spectroscopic and SED classifications in our 
sample is due to the fact that several sources optically classified as normal or starburst 
galaxies do show a Seyfert 2/Seyfert 1.8-like SED, or possibly a composite (starburst+AGN) 
SED. Therefore, the total fraction of type 2 AGN increases from 15.5% (spectroscopic) to 36% 
(SED-fitting), while that of galaxies powered by star-formation (normal+starburst) decreases 
from 70% (spectroscopic) to 48% (SED-fitting). In addition, we have 4% of composite 
(AGN-hstarburst) SEDs that were not considered as spectroscopic class. 

As a partial explanation of these discrepancies, we note that many of the sources spec- 
troscopically classified as normal galaxies do not have enough lines in their optical spectra 
to be properly classified according to the standard line ratios diagnostics (i.e. Tressc et 
al. 1996). In fact, many of them have only Ha and/or [OH] in emission and therefore 
the spectroscopic classification might be misleading. However, in all the cases where the 
classification changed from spectroscopy to SED-fitting, the spectra do either show no AGN 
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emission lines (e.g. [OIII]) at all, or lower than those due to star-formation, making them 
clearly classified as galaxy (normal or starburst) through the diagnostics diagrams. In any 
case, the SED-fitting procedure seems to find more AGN than the optical spectral lines 
classification. One possible reason is that in most of these objects the AGN dominates the 
energetic output only in the MIR. Indeed, most of the sources in our sample with galaxy 
spectroscopic classification and SED reproduced by a Seyfert 2/1.8 template could be fitted 
by a normal galaxy in the optical/NIR part of their spectra, but not in the MIR (in particular 
in the 3-8 fim range), where the data are too fiat to be reproduced by a normal /starburst 
galaxy SED and therefore require the presence of a low-luminosity AGN. The AGN shows up 
just in the range where the host galaxy SED has a minimum, due to the junction between 
the stellar photospheric emission (dominating the optical/NIR) and the dust component 
(starting dominating at A > 5/im; see e.g. the late-type spirals and starburst templates in 
Fig. d]), which is also the range where the hot dust heated by an AGN is expected to start 
increasing towards the peak. In case an AGN component is present with a NIR luminosity 
similar to that of the host galaxy, emission from hot dust heated by the AGN contributes 
to the NIR/MIR filling up the dip observed in star-forming galaxies and producing a fiat 
NIR/MIR spectrum (see e.g. the Seyfert templates in Figured]). This is consistent with 
the results about X-ray selected obscured AGN, which are known to be "elusive" in their 
optical spectra due to host galaxy dilution or heavy dust and gas obscuration or both (i.e. 
see Fiore et al. 2003; Maiolino et al. 2006; Tajer et al. 2007; Caccianiga et al. 2007; Civano 
et al. 2007; Cocchia et al. 2007). Unfortunately, a check with the X-ray luminosity for our 
"elusive" AGN is inconclusive for our purposes, since the XMM observations cover only a 
small portion of the field (15%). 

In order to visualise what discussed above and look at the average properties (in terms 
of redshift and rest-frame luminosity) of the different SED classes, in Fig. H] we show the 
redshift (in logarithmic scale) and the log(Li5^m) distributions for the different SED classes 
to which our sources belong. We notice a clear trend of increasing 15-/im luminosity (and z) 
from early-type (SO, Sa) to later-type (Sd, starburst) SEDs. Type-2 AGN are spread over 
almost the entire range of Lis^m (~10^ - 10^^ L0, except two z > 2 objects at Lis^^ > 10^'^ 
Lq) and over a large redshift range (0.05< z <3), while type I's occupy the region of higher 
luminosities (^15^^ > 10^° L©) and redshifts (most of them are at 0.5 < z < 3). The different 
spectroscopic classifications have been highlighted by filling the relative distributions with 
different patterns, as shown in the legend within the plot. 

Although very uncertain, since the templates often do not fit well the 70 and 160 fim 
data, the total IR luminosities (Ljr) have been computed by integrating the best-fitting 
template SEDs between 8 and 1000 fim. Most of the starburst galaxies and some of the type 
2 AGN are in the ULIG luminosity range (L/r > 10^^ L0), with the remainder in these two 
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classes in the LIG range (10^^ < Lm < 10^^ Lq). The majority of normal galaxies have 
Ljji < 10^^ L0, with just some in the LIG range. All the type 1 AGN are in the ULIG range, 
with few even in the Hyper-LIG (HyLIG) range {Ljji > 10^^ L©), where all the composite 
AGN are indeed. 



5. Discussion and Conclusions 

We have derived the broad-band (from FUV to FIR) SEDs for the largest available highly 

complete (72%) sample of intermediate- 2; MIR selected galaxies and AGN with spectroscopic 
identification and redshift. The majority of these objects are at 0.1 < ^ < 1.3 and in 
the flux density range where strong evolution is observed in the counts of MIR sources. 
Their broad-band SED characterisation is therefore fundamental to understand the nature 
of sources responsible for the observed evolution. Based on the SED-fitting technique we 
have classified the MIR sources, identifying AGN signatures in about 53% of them. This 
fraction is significantly higher than that derived from optical spectroscopy (~29%) and is due 
principally to the identification of type 2 AGN activity in objects spectroscopically classified 
as galaxies. This might be caused, at least in part, by host galaxy dilution in the optical, 
similarly to what is observed for obscured AGN selected in X-ray. It is likely that in most 
of our objects where the AGN luminosity in the optical band is fainter than that of the host 
galaxy, the AGN dominates the energetic output only in the MIR, showing up just in the 
range where the host galaxy SED has a minimum, while the hot dust heated by the AGN is 
expected to start peaking. 

We note that the fraction of MIR sources with an "elusive" AGN derived in this work 
is in agreement with recent results from MIR spectroscopy (with Spitzer-IRS) of local star- 
forming galaxies from the SpitzerSINGS sample by Smith et al. (2007), who find that a high 
fraction (~50%) of local galaxies do harbour low- luminosity AGN (LINER or Seyfert types). 
These low-luminosity AGN, contrarily to the luminous ones which completely destroy PAH 
grains, are able to excite the PAHs and to modify the emitting grain distribution, thus 
providing unusual PAH spectra with very weak or absent 6.2 /im, 7.7 /im and 8.6 /im bands, 
suppressed in strength by up to a factor of 10. An indication of the presence of such weak 
AGN systems, which are found to shift power from the short to the long wavelength PAHs, 
can be revealed, for example, by the strength ratio between the PAH features at 7.7 fim and 
at 11.3 nm (Smith et al. 2007). Therefore, the definitive confirmation of the presence of 
weak AGN systems in our intermediate- 2; MIR dusty galaxies, at present, could be obtained 
only through IRS spectroscopy, allowing to test with fine ratios their SED shape indications. 

Given the new percentages of AGN, we have updated the relative fractions of AGN 
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and galaxies contributing to the observed 15-/im source counts. Although the results of 
our work might still be affected, at some level, by photometric errors in the data and by 
some degrees of degeneracy in the template SEDs, the analysis of the photometric errors 
and the agreement between the primary and secondary solutions in the vast majority of 
the cases make us confident that the majority of the SED-classified AGN are reliable. For 
coherence with the previous classification and source counts (see La Franca et al. 2004 and 
considerations therein) and to avoid biasing towards AGN, we have treated all the sources 
without spectroscopic redshift as galaxies. In fact, a spectroscopic sample 72% complete can 
be biased towards AGN since the large majority of the spectroscopically unidentified objects 
have 15-/im flux densities fainter than 1 mJy. Galaxies dominate over AGN in this range 
of flux and it is therefore likely that most of the unidentified objects are galaxies. On the 
other hand, since we know that optical spectroscopy can miss AGN signatures, we can take 
the fraction of spectroscopic AGN (used for deriving the AGN LF at 15 /xm by Matute et 
al. 2006) as a lower limit. In Fig. [5] we show the fractions of AGN as function of 15-/xm 
flux obtained by considering the spectroscopic classiflcation (dashed line; from Matute et al. 
2006) and the SED classification (filled circles with error-bars; this work). Other recent works 
(Brand et al. 2006; Treister et al. 2006), though based on different diagnostics (e.g. the 24 
fim to 8 yum flux ratio or X-ray detection), estimate lower limits to the AGN contribution in 
MIR surveys signiflcantly higher than those predicted by Matute et al. (2006) and in very 
good agreement with the fractions derived in this paper (see Fig. [5l where the Brand et 
al. 2006 and Treister et al. 2006 results are reported for comparison, with the 24-/im flux 
densities converted to 15 fim by using the opportune 24/15 fim ratio as function of flux as 
computed by Gruppioni et al. 2005). Such results might have a signiflcant impact on the 
interpretation of galaxy and AGN evolution and on the physics of the MIR selected objects, 
so far often erroneously assumed to be all starburst galaxies in many models. In fact, all 
the existing models of galaxy and AGN evolution in the IR should be revised by taking into 
account the higher fraction of AGN dominating in the MIR wave-range. 

By comparing the AGN fraction derived by our SED-fltting analysis to the "lower-limit" 
coming from the spectroscopic classification, we can determine how the relative source counts 
of galaxies and AGN (and consequently the relative evolutionary models) should change 
according to the new values. To this purpose, in Fig. we have plotted the extragalactic 
source counts at 15 fim in the SI field (Gruppioni et al. 2002) with the relative contribution of 
galaxies and AGN computed according to the previous (left) and new (right) determination 
of AGN fractions. The difference between the two is small at bright fluxes, but the SED 
fltting method flnds more and more AGN going to fainter fluxes. The galaxy population 
now dominates below ~3 mJy (previously 6-7 mJy), although now the AGN counts are only 
a factor of 2 lower even at the lower flux densities. For comparison, we have also shown the 
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results of a recent work by La Franca et al. (2007), who updated the AGN source counts at 
15 fim on the basis of X-ray band {XMM) observations (Puccetti et al. 2006) on the central 
0.6 deg^ of ELAIS-Sl. By classifying as AGN all the MIR sources with an unabsorbed 2-10 
keV X-ray luminosity > 10^^ erg s~^, the authors find that at least 13% of the previously 
classified galaxies on optical basis (by La Franca et al. 2004) do harbour an AGN, therefore 
concluding that > 24% of the 15-/im sources with flux density > 0.6 mJy are AGN. The AGN 
source counts at 15 /im updated by La Franca et al. (2007) and plotted in Fig. [6] as dot- 
dot-dashed line are indeed significantly higher that those based on optical classification, and 
in better agreement with (although somewhat lower than) the AGN counts derived in this 
work. The higher efficiency of the SED-fitting technique with respect to X-ray luminosities 
in detecting AGN activity, in this specific case, might be partially due to the fact that XMM 
observations in SI are not deep enough (S(2-10keV) = 2 x 10~^^ erg cm~^ s~^) to detect 
(and thus allow conclusive results for) all (most of) our sources. Moreover, it is known 
that even the deepest X-ray surveys fail in detecting highly obscured (i.e. Compton thick: 
Nh > 10^^cm~^) AGN, while sensitive measurements in the IR range, and in particular SED 
studies in the NIR/MIR range, provide a unique and efficient opportunity to recover the 
fraction of obscured or "elusive" AGN not identified in X-ray surveys (Alonso-Herrero et al. 
2006; Fiore et al. 2008; Martinez-Sansigre et al. 2007; Daddi et al. 2007). 

We have compared our results to those of Hickox et al. (2007), though the sample 
selected by those authors is an IRAC selected one fulfilling also additional selection criteria 
(i.e. detection in all the /i?y4C bands and in the R band, IRAC colours falling in the Stern et 
al. 2005 AGN selection region, spec- or photo-z >0.7), finding that according to the Hickox 
et al. criteria just a few of our AGN can be considered "obscured". In fact, very few of 
our AGN do fall above the R — [4.5] =6.1 separation boundary considered by Hickox et al. 
(2007) to divide "obscured" from "unobscured" AGN. All our type 1 and most of type 2 
are below that boundary, while most of the composite AGN are above. This would imply 
that most of our composite AGN are "obscured" and all our type 1 and most of our type 2 
AGN are "unobscured". The cumulative redshift distribution of our z > 0.7 AGN is in good 
agreement with that derived by Hickox et al. (2007), up to 2; = 1.2 — 1.3. Above z = 1.3 
there are just few sources in our sample, due to our relatively bright 15-/im and i?-band 
selection, therefore we can consider our AGN (in particular type 2) sample complete only 
up to z = 1.3. 

The results of this work will be very useful for updating all the models aimed at inter- 
preting the deep infrared survey data and, in particular, for constraining the nature and the 
role of dust-obscured systems in the intermediate/high-redshift Universe. The main changes 
to the actual evolutionary models for IR galaxies and AGN would consist in better defined 
evolutionary source classes (and SEDs) and in updated evolutions for the different classes. 



- 18 - 



In particular, referring to the specific case of the Pozzi et al. (2004) and Matute et al. (2006) 
models, which are based on the same sample considered here, although star-forming galaxies 
are still the dominant class at low flux densities, their evolution should be lower than that 
derived by Pozzi et al. (2004). Similarly, the AGN (mainly type 2's) should evolve more 
rapidly (and probably both in luminosity and density) than derived by Matute et al. (2006), 
similarly to star-forming galaxies. In fact, the higher fraction of AGN in infrared surveys 
than found through optical spectroscopy is a result common to different works (i.e. Brand 
et al. 2006; Treister et al. 2006) and applies mainly to type 2 AGN. Type 1 AGN are in fact 
easier to identify in optical, due to their broad lines, thus most (all) of them are classified 
from their optical spectra. Type 2 AGN are more "elusive" and difficult to reveal in the 
optical. Therefore, the higher fraction of IR AGN is mainly due to the unveiling of type 2 
AGN previously erroneously classified as galaxies. For this reason, the starburst galaxies, 
to whom all the evolution observed in the MIR was commonly attributed, are probably in 
a smaller number than formerly believed and part of their evolution should be ascribed to 
type 2 AGN. This is also evident from the right panel of Fig. [6l where the revised AGN 
source counts increase towards fainter fiuxes with a slope similar to (though slightly lower 
than) that of galaxies, while, according with the results based on optical spectroscopy (left 
panel), the AGN were rapidly converging towards the lower fiux densities. 

Finally, instead of considering non-evolving normal galaxies (and therefore ascribing all 
the evolution to starburst galaxies), it would be more appropriate to consider galaxy SEDs 
evolving with z (and/or V) from early- type to late-type (and to starburst), as suggested 
by the results shown in Fig. HJ and eventually different evolutions for different luminosity 
intervals. The revision of evolutionary models for IR sources based on the results presented 
here is beyond the aims of this work and will be treated in a forthcoming paper (Gruppioni 
et al., in preparation). 

To summarise, in this paper through the broad-band SED analysis of the largest avail- 
able highly (72%) complete spectroscopic sample of MIR selected galaxies and AGN at 
intermediate z, we have: 

• verified the assumption that local template SEDs are able to reproduce also most of 
the MIR-selected galaxies and AGN at 0.1 < z < 1.3; 

• found AGN activity in a significantly higher fraction of sources (~53%) than derived 
from optical spectroscopy (~29%); 

• derived new relative fractions of AGN and galaxies contributing to the observed MIR 
source counts, with AGN going from ~10-20% of the MIR population at 6*15 <0.6 mJy 
up to ~ 90 - 100% at 5i5 > 10 mJy; 
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computed new source counts in the MIR, with the AGN counts (especially type 2's) 
now increasing with a slope similar to that of galaxies. This result is likely to imply 
different evolutionary rates than formerly considered by the present models (i.e. lower 
rates for starburst galaxies and higher rates for type 2 AGN). 
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Table 1. Multiwavelength Properties of the 15 //m Sources in the ELAIS-Sl Field: Example Table 
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Table 2. Spectroscopic Classification versus SED classification 



SED ^ \ SPEC 


Tot 


AGNl 


AGN2+LINER 


Starburst 


Galaxy 


Unclass 


Tot 


203 


27 


25+7 


32 


107 


5 


AGNl 


24 


20 (74.1%) 


1+0 (3.1%) 


2 (6.2%) 


1 (0.9%) 





AGN2 


72 


3 (11.1%) 


14+4 (56.3%) 


14 (43.8%) 


35 (32.8%) 


2 (40.0%) 


AGN+Starburst 


11 


4 (14.8%) 


5+0 (15.6%) 


1 (3.1%) 


1 (0.9%) 





Starburst 


13 





2+1 (9.4%) 


6 (18.8%) 


3 (2.8%) 


1 (20.0%) 


Galaxy 


83 





3+2 (15.6%) 


9 (28.1%) 


67 (62.6%) 
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Fig. 1. — Library of template SEDs from Pollctta ct al. (2007). The SEDs arc plotted in 
arbitrarily scaled luminosity vSv versus wavelength. 
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Fig. 2. — Example of SED fits for our sample: the data (filled circles with error-bars), plotted 
as flux {Si, in /xJy) versus observed wavelength, are over-imposed to the corresponding best- 
fitting template SEDs (sohd hne), opportunely redshifted at the source's z. The plotted 
uncertainties are those reported in the various catalogues and not those used for the fits. 
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Fig. 3. — Rest-frame SEDs (dots) of 15-fim sources belonging to the different broad SED 
classes (AGNl, ULIG, AGN2, Starburst and Spiral galaxy) compared to the most frequent 
best-fit template SED of that class, all normalised to the iCs-band flux density. The white 
triangles represent the average values that the templates should have in the different rest- 
frame bands to better reproduce our sources, while the grey-shaded areas show their relative 
uncertainty region corresponding to ilcr. 
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Fig. 4. — Logarithmic distributions of z (left) and Lis^^ luminosity (right) for the different 
SED classes reproducing our data. The luminosities have been derived by considering for 
each source its best-fitting template SED. Within each distribution the different spectroscopic 
classes are highlighted, as explained in the legend. 
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Fig. 5. — Fraction of AGN over total extragalactic sources as a function of the 15-/xm flux 
density, as derived from the spectroscopic classification (dashed line; Matutc ct al. 2006) and 
from the SED-fitting performed in this work (filled circles with error-bars). For comparison, 
the results of Brand et al. (2006) and Treister et al. (2006) have been plotted as open 
squares and triangles (open: non corrected for AGN not detected in X-rays; filled: corrected 
for AGN not detected in X-rays), respectively. 
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Fig. 6. — Observed extragalactic source counts at 15 /im in ELAIS-Sl: total counts with 
uncertainties (dark-grey shaded area); type 2 AGN contribution (hght-grey shaded area); 
type 1 AGN contribution (dotted area); galaxy contribution (hatched area). The total (type 
1 + type 2) AGN contribution is shown by the dashed line, while the dot-dot-dashed line 
represents the total AGN counts at 15 /xm updated by La Franca et al. (2007) on the basis of 
hard X-ray observations. In the left panel the AGN/galaxy contributions are derived from 
the spectroscopic classification (La Franca et al. 2004), while in the right panel they come 
from the SED-fitting work described in this paper. All the optically unidentified sources 
have been conservatively considered as galaxies in both plots (see Section [S]). 



